MZ-TH/01-06 



hep-ph/0103047 
March 2001 



Low-energy gluon contributions to the 
vacuum polarization of heavy quarks 

S. Groote^ and A. A. Pivovarov^'^ 

^ Institut fiir Physik der Johannes-Gutenberg-Universitat, 
Staudinger Weg 7, 55099 Mainz, Germany 

^ Institute for Nuclear Research of the 
Russian Academy of Sciences, Moscow 117312 



Abstract 

We calculate a correction to the electromagnetic current induced by a heavy quark loop. 
The contribution of this correction to the vacuum polarization function appears at the 
O(a^) order of perturbation theory and has a qualitatively new feature - its absorptive 
part starts at zero energy in contrast to other contributions where the absorptive parts 
start at the two-particle threshold. Our result imposes a constraint on the order n of the 
moments used in the heavy-quark sum rules, n < 4. 
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High precision tests of the standard model remain one of the main topics of particle phe- 
nomenology. The recent observation of a possible signal from the Higgs may complete the 
experimentally confirmed list of the standard model particles Because experimental data 
are becoming more and more accurate, the determination of numerical values of the parameters 
of the standard model lagrangian will require more accurate theoretical formulae. Recently an 
essential development in high-order perturbation theory calculations has been observed. A re- 
markable progress has been made in the heavy quark physics where a number of new physical 
effects have been described theoretically with high precision. The cross section of top-antitop 
production near the threshold has been calculated at the next-to-next-to-leading order of an 
expansion in the strong coupling constant and velocity of a heavy quark with an exact account 
for Coulomb interaction (as a review, see Ref. [Q). This allows for the best determination of 
a numerical value of the top quark mass. The method of Coulomb resummation resides on a 
nonrelativistic approximation for the Green function of the quark-antiquark system near the 
threshold and has been successfully used for the heavy quark mass determination within the 
sum rules ^, Being applied to quarkonium systems this method is considered to give 
the best estimates of heavy quark mass parameters P, 1^, ^ . Technically an enhancement of 
near-threshold contributions to sum rules is achieved by considering integrals of the spectral 
density of the heavy quark production with weight functions which suppress the high-energy 
tail of the spectrum. The integrals with weight functions for different positive integer n, 
s = E"^, where E is the total energy of the quark-antiquark system, are called moments of the 
spectral density and most often used in the sum rules analysis [|10 . 



In the present note we show that there is a strong constraint on the order n of the moment 
that can be used in heavy quark sum rules. Because of the contribution of low-energy gluons, 
only moments with n < 4 exist if theoretical expressions for the corrrelators include the 0{al) 
order of perturbation theory. 

The basic quantity for the analysis within sum rules is a vacuum polarization function 11 (g^) 

Un'z J {Tj,{x)UO))e''''d'x = {q,q, - g,.q'Mq') (1) 

of the vector current = q'-f^q of a heavy fermion q. With the spectral density p{s) defined 
by the relation 

p{s) = ^(U{s + iO)-U{s-iO)), s>0 (2) 
Zm 



the dispersion representation 

,2 



n(,") = / ^ (3) 

s — q"^ 



holds. A necessary regularization (subtractions, for instance) is assumed in Eq. (|^). The 
normalization of the vacuum polarization function Il{q'^) in Eq. (|I|) is chosen such that one 
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obtains a high-energy hmit hm^^oo p{s) = 1 for a lepton. For the quark in the fundamental 
representation of the SU{Nc) gauge group a high-energy hmit of the spectral density reads 
p(oo) = Nc- The integral in Eq. runs over the whole spectrum of the correlator in Eq. (|l|) 
or over the whole support of the spectral density p(s) in Eq. (Q). The moments of the spectral 
density p{s) of the form 

(4) 

are usually studied within the sum rules method for heavy quarks [0. These moments are 
related to the derivatives of the vacuum polarization function 11 (g^) at the origin, 



(5) 



Such moments are chosen in order to suppress a high energy part of the spectral density p{s) 
which is not measured accurately in the experiment. Within the sum rule method one assumes 
that the moments in Eq. @ can be calculated for any n or, equivalently, that the derivatives in 
Eq. @ exist for any n. The existence of moments seems to be obvious because one implicitly 
assumes that the spectral density of the heavy quark electromagnetic currents p{s) vanishes 
below the two-particle threshold s = 4m^, which means that the vacuum polarization function 
of heavy quarks n(g^) is analytic in the whole complex plane of except for the cut along 
the positive real axis starting from 4m^. This assumption about the analytic properties of 
the vacuum polarization function n(g^) is known to be wrong if a resummation of Coulomb 
effects to all orders of perturbation theory is performed: as a result of such a resummation 
the Coulomb bound states appear below the perturbation theory threshold s = 4m^. The 
assumption that the moments in Eq. (^) exist for any n may also be wrong in high orders 
of perturbation theory in models with massless particles, for example, in QCD with massless 
gluons. The validity of this assumption depends on details of the interaction. In QCD, at the 
0{al) order of perturbation theory there is a contribution of massless states to the correlator 
in Eq. (|T]) that leads to the infrared (small s) divergence of moments for large n because of the 
branching point (cut) singularity of n(g^) at the origin. We determine the behaviour of the 
vacuum polarization function n(g^) at small (g^ <^ m^) as 

n(g^)|,..o = C^.(^)'ln(^) (6) 
with ^ 

Here date are the totally symmetric structure constants of the SU{Nc) gauge group defined 
by the relation dabc = 2tr({t", and are generators of the group with normalization 
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Figure 1: Heavy quark loop correction to the electromagnetic current 



tr(t"t^) = 1/2. For the SU{3) gauge group of QCD one has dabcdabc = 40/3. The parameter 
in Eq. (|^) is the renormalization point. 

The singularity of the vacuum polarization function given in Eq. (a cut along the positive 
real axis in a complex g^-plane) prevents one from calculating moments of the spectral density 
in Eq. (§) with n > 4. Indeed, the high order derivatives of n(g^) at the origin determining the 
high order moments according to Eq. do not exist for n > 4 because of a branching point 
singularity as one can see from Eq. (^. In terms of the moments one can see this by calculating 
the behaviour of the spectral density at small energy s, 

p{s)U.o = C,[£,y (8) 

which makes integrals in Eq. divergent at small s for n > 4. The formulae for the vacuum 
polarization function in Eqs. @ and (0) are given for a heavy quark in the SU{Nc) (E) U{1) 
gauge model. The result for QED may be obtained with the obvious substitution a for 

the coupling constant and by setting dabcdabc = 1- Contributions of light (massless) quarks 
appear in the 0(af) order of perturbation theory and are neglected. 

We now present a derivation of our result given in Eqs. (^) and (0) and discuss some 
consequences for the phenomenology of heavy quarks. Note that the induced current is a 
correction of order 1/m^ in the inverse heavy quark mass which vanishes in the limit of an 
infinitely heavy quark. Corrections in inverse heavy quark masses are important for tests of 
the standard model at the present level of precision and have been already discussed in various 



areas of particle phenomenology [0, |T^, |T3|. 

A diagram which gives a correction to the electromagnetic current due to a heavy quark 
loop is given in Fig. [l|. Two-gluon transitions are forbidden according to a generalization of 
Farry's theorem to nonabelian theories |T^. We are interested in a behaviour of the amplitude 
associated with the diagram Fig. |I] at low energy and therefore take the limit of a very heavy 
quark. Formally the limit m ^ cxd is taken which in physical terms means that m is much larger 
than all momenta of external legs of the diagram, namely the three gluons and the photon. The 
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induced current is written in a covariant form as a derivative of an antisymmetric operator 
C^jy built from the gluon fields only, 

= d,0^,, O^, + = 0. (9) 

This structure of the induced current automatically guarantees the current conservation 

d^J^ = (10) 

as it should be for the electromagnetic current. 

A straightforward calculation of the diagram presented in Fig. ^ gives the result for the 
induced correction to the electromagnetic current 

^'^ = Ii4^ (5 W + 145.0n (11) 

with 

fnlJ-i^ ^ ria r^h r^c /^A"^ A n°- H'^ 

where is a gauge field strength tensor for the gauge group SU {N^) . 
A correlator of the induced current J'^ has a general form 

{TJ^{x)J,{0)) = -d^dp{TO,MOup) (13) 

where an explicit expression of the current as a derivative of the antisymmetric operator O^u 
has been employed. The resulting correlator {TOfj_a{x)Oui3) in Eq. ([T3|) contains only gluonic 
operators as is seen from Eqs. (0) and (|12]). Such correlators were considered previously in 
the framework of perturbation theory |T5| , |TB[ . In the leading order of perturbation theory the 
correlator in Eq. ( |T3|) has a topological structure of a sunset diagram. Technically, a convenient 



procedure of computing the sunset-type diagrams is to work in the configuration space [17 
We find 

{TJ,{x)U^)) = -^3^^ - g,.d-) -L (14) 

A Fourier transform of the correlator in Eq. ([1^) gives the vacuum polarization function in 
momentum space which reads 

U-n^ I {TJ,ix)MO))e''''d'x = Cg{q,q, - g,.q') In (^-^^ (15) 

with the constant Cg taken from Eq. (|^). The spectral density of the polarization function in 
Eq. (|1^) is given in Eq. (||). 

Note that the spectral density of the correlator in Eq. (|I^ can be found without an explicit 
calculation of its Fourier transform. Instead one can use a spectral decomposition (dispersion 
representation) in configuration space. 



with D(a;^, s) being the propagator of a scalar particle of mass ^/s, 



2 2. _ %m^r^K^{m4^) 

where Ki{z) is a McDonald function (a modified Bessel function of the third kind, see e.g. 
Ref. [1^). r(2) is Euler's gamma function. 

An asymptotic behaviour of the spectral density of the corresponding contribution for large 
energy (when the limit of massless quarks can be used) is well known [0, [2^, |2l| and reads 

'as\^ dabcdabc / 176 ^, , 

nj 1024 V 3 ' 

Here ({z) is the Riemann ( function and C(3) = 1.20206 . . .. The contribution to the spectral 

density given in Eq. ([I8|) is negative while our result given in Eq. (H) is positive as it should be 

for the spectral density of the electromagnetic current which is an Hermitian operator. 

In QCD we find 

Our result has an immediate application to the determination of heavy quark parameters 
within the method of sum rules. Because of the low-energy gluon contributions, the large n 
{n > 4) moments of the spectral density do not exist and cannot be used for phenomenological 
analyses. Note that in early considerations of sum rules quite large n were used. For instance, 
the numerical value of the gluon condensate was extracted from sum rules for the moments with 



n ~ 10 -i- 20 |110| , P^ . In view of our result one has either to limit the accuracy of theoretical 
calculations for the moments to the 0(0;^) order of perturbation theory which seems insufficient 
for a high precision analysis of quarkonium systems (especially if the Coulomb resummation in 
all orders is performed) or to use only a few first moments with n < 4. For small n, however, 
the high-energy contribution, which is not known experimentally with a reasonable precision, 
is not sufficiently suppressed and introduces a large quantitative uncertainty into sum rules 
for the moments. An alternative analysis based on finite energy sum rules is free from such a 
problem and can be used in phenomenological applications 



Note in passing that there is no low-energy gluon contribution (and low-energy divergence 
problem) for correlators of the currents containing only one heavy quark with mass m. The spec- 
trum of such correlators starts at and there are no massless intermediate states contributing 
to the correlator in perturbation theory. The theoretical expressions for such correlators can be 
used for high precision tests when the accuracy of experimental data in correspondent channels 
will improve in the future. 

To conclude, we have presented a correction to the electromagnetic current of a heavy quark 
induced by a virtual heavy quark loop. The spectrum of the correlator of such an induced 
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current starts at zero energy. This fact makes impossible the standard analysis of the moment 
sum rules for n > 4 at the 0{a^) order of perturbation theory. 
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